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Abstract 

■ In this paper we present a self-contained variational theory of the layer potentials 

for the Stokes problem on Lipschitz boundaries. We use these weak definitions to 

show how to prove the main theorems about the associated Calderon projector. 

Finally, we relate these variational definitions to the integral forms. Instead of 

working these relations from scratch, we show some formulas parametrizing the 

Stokes layer potentials in terms of those for the Lame and Laplace operators. While 

Q\ ■ all the results in this paper are well known for smooth domains, and most might be 

known for non-smooth domains, the approach is novel a gives a solid structure to 

^ ' the theory of Stokes layer potentials. 
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1 Motivation and basic notation 

In this paper we present a coherent and self-contained theory of the layer potentials for 
the Stokes equation on Lipschitz boundaries, in two and three dimensions. The core 
of the paper is the presentation of the results in variational form, using a formalization 
based on weighted Sobolev spaces. This presentation, a la Nedelec, follows the way how 
the properties of integral operators/equations (in particular coercivity of equations of 
the first kind) were introduced in the often quoted and not easy to find lecture notes of 
Jean-Claude Nedelec [13] . What is missing there is a rigorous proof that the potentials 
correspond to the integral expressions that have been known of old (with some work, this 
recognition can be done for smooth enough domains). A way to solve this problem is to 
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start from the other side, defining the potentials in integral form, and then showing that 
they satisfy the variational equations that were used for the weak definition of potentials. 
The definition of potentials, in weak form but using integral expressions, is the modern 
standard for their introduction, or better said, it is the standard that numerical analysts 
have accepted as the one that satisfies their theoretical needs. A big part of the boundary 
integral community has embraced William McLean's monograph [13] as the answer to their 
prayers, the book that proves everything, the source of all needed references. However 
general it is, the general frame that is developed in [T3] does not include the Stokes 
equation. We are actually going to openly borrow results from this text, in a slightly 
surprising way, at the very end of this paper. The long awaited -and finally appeared 
in the inmensity of its six hundred page format- book of George Hsiao and Wolfgang 
Wendland [TT] deals with many more equations, but in its desire to enjoy the full power 
of the pseudodifferential calculus, it does not include many results (and/or tools) to deal 
with layer potentials on, say, polyhedra. This fast review would be definitely incomplete 
if it were not to mention a third way, which has proven to be the one providing the best 
mapping properties for boundary integral operators. The techniques employed in this 
way are much more sophisticated than what we are going to see in this paper, picking 
from advanced analysis of singular integrals. The Chicago-Minnesota school of harmonic 
analysis (Eugene Fabes, Carlos Kenig, Gregory Verchota [5]) started the development of 
this kind of treatment of the Stokes integral operators. 

This is what we are going to do next. We are going to use the variational technique 
of Nedelec to introduce the layer potentials on Lipschitz surfaces for the Stokes problem 
in two and three dimensions and we are going to use these definitions to prove the main 
properties of the associated boundary integral operators. Next, we are going to recon- 
struct the potentials using integral formulas. With combinations of Lame and Laplace 
potentials, plus the standard properties thereof [13j . we will show that integral potentials 
and variational potentials are one and the same. 

Some (or even all) of the results of this paper are known, but not well known, and 
definitely not easy to find in the literature. The presentation is original (up to the fact that 
we mimick the French school variational method) and we make a point in solving some 
difficulties in the two dimensional case that otherwise lead to a definition of potentials 
in a quotient space. To the best of our knowledge, we are novel in integrating the two 
approaches in one unique proof format. In any case, the goal of the paper is twofold: give 
a clean presentation of techniques and easy-to-refer results; help the reader get acquainted 
with the mathematics behind the theory of layer potentials. 

Apart from well-established background in basic functional analysis, Sobolev space 
theory and elliptic PDE, we will need some popular results (Korn's inequalities and the 
theory of abstract mixed problems, a.k.a., the Brezzi conditions) and a few more ingredi- 
ents: 

• Of the big collections of weighted Sobolev spaces developed by Bernard Hanouzet 
and Jean Giroire, we will employ one of each, with the aim of controlling behavior 
at infinity of potentials with local H 1 regularity. The results that are needed are 
taken from [21 HP]. 

• A result concerning the surjectivity of the divergence operator from the previous 
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weighted Sobolev spaces to L 2 will be used. This was proved by Vivette Girault 
and Adelia Sequeira [9]. 

• Finally, the integral theory of layer potentials for the Laplace and Lame operators 
will be used at some key points at the end of this paper. The results will be taken 
from William McLean's monograph [13], which exploits the ideas of a celebrated pa- 
per by Martin Costabel [5]. (Note again that the Stokes operator is not a particular 
case of the very general theory developed in [13"].) 

Precise reference on all the results will be given as we reach them. We next list the basic 
background material and the general notation that will be used throughout. All elements 
of this list are standard and the reader can use this collection of bullet points for easy 
reference. 

Functional analysis. If if is a Hilbert space, H := H d endowed with the product norm: 
d £ {2, 3} will be space dimension. If if is a Hilbert space and X C H, then H' is 
the dual space of H and X° = {£ £ H' : £(x) = Wx £ X} is the polar set of X. 

Scalars, vectors, and matrices. We will keep separate notation for scalar (italic, as in 
u), vector (boldface, as in u) and matrices-tensors (capitals, as in U) quantities. 
In case of doubt, vectors are column vectors. The colon will be used for the tensor 
(Frobenius inner) product of matrices U : V = X/yUyVy. Everything will be 
real-valued and there will be no need for conjugation. The tensor product of two 
non-zero vectors u and v is the rank one matrix u <g> v with entries U\ Vj and the 
following formula should be kept in mind: (u <g> v)w = (v • w) u. 

L 2 products. Parentheses will be used for L 2 products: 

(u,v) c := [ uv, (u,v) = / uv, (U,V) C := / U : V. 
Jo Jo Jo 

The associated norms will be denoted || • \\o- 
Multivariable calculus. Given a vector field u, we define Du and e(u) by 

i I ><«),,: j£, 6(u):=i(Du+(Du) T ). 

Note the useful formula e(u) : e(v) = e(u) : Dv. If A is matrix valued, then divA 
is the divergence operator applied to the rows of A (so that the output is a vector). 
This is consistent with the formula 

div(A T v) = (div A) • v + A : Dv, 

which is in the heart of the integration by parts process. 

The geometric setting. In what follows f2_ is a connected bounded Lipschitz open set 
boundary T. The exterior domain is Q + := M d \ The case of boundaries with 
multiple connected components is very easy to handle, but we will prefer to stay in 
the connected case. 



3 



Sobolev theory. For a bounded or unbounded domain on one side of its Lipschitz bound- 
ary, H 1 ^) is the usual Sobolev space, 7 : — > H l l 2 (dO) is the trace operator, 
and Hq(0) = Ker7. When we identify L 2 (dO) with its dual space, the dual space of 
H l l 2 (dO) is automatically represented by a superset of L 2 (dO). This representation 
is called E~ x l 2 {dO). The duality product H- l / 2 (dO) x H l ' 2 (dO) is denoted (•, -) do . 
For matrix valued A G H l l 2 (dO) dxd and £ G H~ l / 2 (dO), we will use (A,£) ao to 
denote the vector with components ^j=i(£j' ^ij)do- The space of C°° functions with 
compact support contained in O will be denoted V(0). 

Local behavior. The subscript loc will be used in some parts of this document. For 
instance, when we write u G H^ oc (0), we mean that ipu G H l {0) for all 99 G 
(not only </? G X>((9)). This means that we are only localizing at infinity and we are 
not separating ourselves from the boundary of the set. 

2 The Sobolev spaces of Giroire and Hanouzet 

For an open set O C M. d , we consider the weighted Sobolev space 

W(0) := {u : O -> R : pu G L 2 (C), Vm G L 2 (C)}, 
endowed with its natural norm 

\\ u \\l,p,o '■= \\P u \\o + l|V«||o, 

where p = pa is given by: 

, , 1 r , 1 1 

V 1 + l x l 1 + 2 ^Sl 1 + l x l ) V 1 + l x l 

Note that if O is bounded H l (0) = W(0). It is an easy exercise to prove that if O is 
unbounded and Pk(0) is the space of rf-variate polynomials of degree no larger than k, 
then 

that is, W{0) does not contain non-trivial polynomials in three dimensions, but it contains 
constant functions in two dimensions. 

Proposition 2.1. The space T>(Jxi d ) is dense in W(M. d ). Moreover, the following norms 
are equivalent to the norm || ■ \\i iPi Rd in W{R d ): 

(a) || V • || R s (for d = 3). 

(b) (||V-|| 2 2 + |j(-)| 2 ) 1/2 (ford = 2). 

In (b), j : W(M. 2 ) — > R is a bounded functional such that j(l) 7^ 0. 
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Proof. The density result in the case d = 3 is easy [TO] Theoreme 1.1]. The case d = 2 
can be found in pQ: it requires to first show that H l (M. 2 ) is dense in ^(R 2 ) using a clever 
damping argument. Part (a) is proved in [TU| Theoreme 1.2]. Part (b) can be found in [TJ, 
stated in equivalent form in a quotient space: ||V- ||r2 is an equivalent norm in W(R 2 )/Po. 
Note that the references for the two dimensional results are more recent than the results 
themselves, which were already mentioned in [12]. □ 



Corollary 2.2. The following norms are equivalent to the product norm in W(R d ): 

(a) ||£(-)I|k 3 (ford = 3). 

(b) (IK")II^ + W-)| 2 ) 1/2 (ford = 2). 

In (b) j(u) := (j(ui) , j(u2)) and j : iy(R 2 ) — > R is a bounded functional such that ^ 0. 
Proof. For all u G T>(R d ) we can write 

2(e(u), e(u)) R d = 2(e(u), Du) K d = (Du, Du) Rd + (divu, divu) Rd . 

Therefore 

||Du|| 2 d < 2||£(u)|| 2 d = |||Du+ (Du) T || 2 d < 2||Du|| 2 d Vu G T>(R d ). 
By density (Proposition 12. ip it follows that 

^||Du|| R d < ||e(u)|| R d < ||Du|| R d Vu e W(R d ). (2.2) 

The result is then a direct consequence of f)2.2p and Proposition 12.11 □ 

Because the spaces W(0) are locally if 1 ((9), we can define the interior and exterior 
traces 

7+ : W(n + ) -> H 1/2 (T), 7 " : H\n_) -> H l ' 2 {T). 
These trace operators are surjective. 

Proposition 2.3. The space T>(Q + ) is dense in Wq(Q+) := {u G W(Q + ) : 7 + m = 0}. 

Proof. This result is a straightforward consequence of Proposition 12.11 and the fact that 
T>{0) is dense in the kernel of the trace operator for any Lipschitz domain O. □ 

Proposition 2.4. The divergence operator div : W(IR d ) — > L 2 (M. d ) is surjective. 

Proof. The divergence operator div : Wo(fi+) — > L 2 (Q + ) is surjective. This is [9] Theorem 
3.3]. Given p G L 2 (M. d ) we takep_ := p\n_ and choose u G H 1 (fi_) such that divu_ = p_ 
in Let ul G H x (R d ) be such that u!\q_ — u_. We then build v + G W (fi+) such 
that divv + = p\n — divul|^ + . Finally, we consider 



u : = 



u_, in fl_ 

ul| n+ +v + , in fi H 



Then u G W(R d ) and div u = p. □ 
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Consider the spaces Ai := Aid'- 

M 3 := {a + b x x : a, b G R 3 }, 

M 2 := {a + 6 (x 2 , —Xi) : a G R 2 , b G R}. 

Note that if v G Ai, then e(v) = and divv = 0. Also dimA^ 3 = 6, dim Ai 2 = 3, and 
Ai fl HQ(fi_) = {0}. Therefore, the space Afr := lAi has the same dimension as Ai and 
elements of Ai are determined by their values on T, which will allow us to just identify 
Air = Ai. Finally, by flZ2) 

Mnw w = (nw, h = 2! < 2 - 3 > 

Lemma 2.5 (Rigid motions and Korn's inequality). Lei O be a bounded domain. 

(a) 7/v G H 1 ((9) satisfies e(v) = an<i (9 zs connected, then vgM. 

(b) There exist constants C, c > snc/i t/iat 

||e(u)|& + c||u|||, > C||Du||^ Vu G H x (0). 

(c) If m := dimA^d and j : H 1 ((9) — >■ R m is linear, bounded, and such that j(c) ^ 
for all 7^ c G Aid, t/ien 

|s(u)||^ + |j(u)| 2V: ' 



an equivalent norm in H 1 ((!?). 

Proof. The proof of Part (a) for d = 3 is given as [131 Lemma 10.5]. The two dimensional 
case is a simple consequence of the three dimensional case, by immersing the two dimen- 
sional vector field in a space of three dimensional vector fields. Part (b) is the well known 
Korn's second inequality [T3], Theorem 10.2]. Part (c) follows from (b) and a compactness 
argument. □ 

Lemma 2.6. If u G W(fi + ) satisfies e(u) = 0, then u G Ai H W(f2 + ). 

Proof. Consider a domain D := £>(0; i? ) \ ^- = -8(0; i?o) H fi+ for large enough, and 
consider the annular domains D n :={x : Rq + n < |x| < Rq + n + 1}. Then, by Lemma 



I2.5f a). there exists c n G AI such that u|o n ee c ra . Comparing the traces of c n _i and c n 
on dB(0; Rq + n), we show that c n ee c n+ i. Therefore u ee c G AI in fi + . □ 

Proposition 2.7. T/ie following norms are equivalent to the product norm in W(f2 + ). 

(a) |H-)lk (ford = 3). 

(b) (lk(')ll^ + + |j(-)| 2 ) 1/2 f/ord = 2j. 

In (b) j(u) := (j(m 1 ), ]{u 2 )) andj : ^(R 2 ) — >■ R zs a bounded functional such that ^ 0. 
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Proof. Take <p G P(R ) such that = 1 in a neighborhood of and </? is supported 
inside a ball 5. Decomposing u = (1 — ip)u + </?u, we can understand (1 — </?)u as an 
element of W(R d ), while ipu is compactly supported in a ball 5. We let B + := B D 

We can apply Corollary 12.21 using j(u) := J r 7^ for the moment being. Simple com- 
putations show then that 

IM|i, ft n+ < 11(1 - ^ulli^Rd + ||yu||i ii0iB+ < Ci||e((l - y)u)|| R d + C 2 ||<y3u||i )jB+ 
< Ci||e(u)||n + + C 3 ||u||b + + C 4 ||u||i )S+ . 

Here we have used that — (p)ui) = and 

£((1 - (p)u) = (1 - <p)e(u) - ~u <g> V<f - ® u 

Applying Korn's inequality (Lemma 12.5T b)) in B + , we can bound 

||u||i,b + < C5^||e(u)||fl + + \\u\\ B j 

and therefore 

IMIwV < C 6 (||e(u)||^ + + ||u||| 4 

The remaining part of the argument is slightly different in two and three dimensions: 

(d = 3) By Proposition 12.61 and (12. 3p . ||e( • )||n + is a norm in W(fi+) which admits the 
usual norm as an upper bound. On the other hand the compact term || • \\b + 
(note that W(B + ) = iJ 1 (5 + ) is compactly embedded in L 2 (B + )) makes this norm 
equivalent to the usual one. This means that we can take this compact term out of 
the norm and we still have an equivalent norm. 

(d = 2) By Proposition EOS and O, (||e( • + L?(-)| 2 ) 1/2 is a norm in W(fi + ) and we 
can bound 

||u|| 1)P ,n + < C 6 (\\e(u)\\l + + |j(u)| 2 + ||u||| + ) 1/2 < C 7 \\u\\ hPj n + . 

Using the same argument as in the three dimensional case, we can eliminate the 
|| • ||_b + term and still get an equivalent norm. 

□ 

3 Review # 1: normal stresses 

Consider a matrix valued function cr (stress) such that 

creLUR d ) dXd , div<reLL(K d \r). 

We can now choose an arbitrary ball B containing f2_ and write B + := Q + (IB. Then we 
define: 

(<r-n,7v) r := (cr, Dv) n _ - (div cr, v) c _ Vv G H 1 (fi_), 
(cr + n,7v)r := -(*, Dv) B+ + (div cr, v) B+ Vv G Hj(S). 



These are matrix extensions of the classical definitions of normal traces in H(div, fi_) 
and H loc (div, With pretty much the same proof as in the scalar case, you show that 
o^n G H _1 / 2 (r) and that we can bound 

||<7- _ n||_i/ 2j r < C(||<r||n_ + ||divcr|| n _), ||cr~n||_ 1/2i r < C(||cr|| B+ + ||div<x|| B+ ). 

We can also write [crn] := cr n — cr + n, so that 

([an],7v) r = (cr, Dv) Rd \ r - (div cr, v) M d\ r Vv G Hj(S). 

If cr G L 2 (R d ) dxd , and diver G L 2 (M d \ T), then we can apply a density argument and 
show that 

([<rn],7v) r = (cr, Dv) Rd \ r - (div cr, v) Rd \ r Vv G H 1 ^). 

Example # 1: Laplace. If u G Hj 0C (R d \ T) and Au = in R d \ T, choosing 

cr A (u):=Du, (3.1) 

we get to the integration by parts formula applied to the components of u 

([ ( 7 A n], 7 v) r = ([<9 n u], 7 v) r = (Du, Dv) K . Vv G Hj(fl). (3.2) 

This same definition can be applied for the non-homogeneous case, when A-tu G L 1 2 oc (lR d \ 
r) (by A± we are meaning the Laplace operator applied in Q ± , not in R d ). 

Example # 2: Lame. If u G Hl c (R d \ T) , we define 

cr L (u) := 2/i£(u) + A(divu)I (3.3) 

and we assume that divcr L (u) = in R d \T (that is, u is a solution of the Lame equations 
in R. d \ T), then we arrive to Betti's formula 

([cr L (u)n], 7 v) r = (cr L (u),e(v)) RAr (3.4) 
= 2//(e(u),e(v)) R d\ r + A(divu, div v) Rd \ r Vv G Hj(B). 

Example # 3: Stokes. If (u,p) G Hj )C (R d \ T) x L 2 oc (M d ), we define 

cr s (u,p) := 2i/e(u) -pi, (3.5) 

and assume that div cr s (u, p) = in R d \ V, we obtain a definition of jump of the normal 
stress for the Stokes problem. Since this is a paper on the exterior Stokes problem, 
let us detail what we obtain in this case. The interior and exterior normal stresses for 
homogeneous solutions of the Stokes problem 

(u,p) G W(M d \T) x L 2 (R d ) such that - 2z/dive(u) + Vp = in R d \ T, 
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are then given by t ± (u,p) = cr s (u,p)n, i.e., 

(t-(u,p), 7 v) r = 2i/(e(u),e(v)) n _-(p,divv) n _ Vv G H 1 (Q_), 
(t + (u,p), 7 v) r = -2^( e (u), £ (v)) n+ + (p,divv) n+ VvGW(Q + ). 

(Note how we have found it convenient to extend the exterior test space to the corre- 
sponding weighted Sobolev space.) We also have 

||t ± (u,p)||_ 1/2 , r < C7 r (||e(u)|| n± + ||p||n ± ). (3.6) 

Proposition 3.1. 

(a) For smooth enough (u,p), t ± (u,p) = 2z/7 ± (e(u)) n — (7 ± p)n. 

(b) For no interior velocity and constant pressure, the normal stress is proportional to 
the normal vector field: t~(0, 1) = — n. 

(c) In dimensions d = 2 and d = 3, t~(u,p) G H 1 ^ 2 (r). 

(d) When d = 2, t + (u,p) G H 1/2 (r). 

Proof. Parts (a) and (b) follow from the definition. For part (c), take v = e G M d . For 
(d), do the same thing in the exterior domain. □ 

4 Thinking of the boundary 

An issue with decompositions of vector fields on a Lipschitz boundary is related to the fact 
that the normal vector field is allowed to have discontinuities and therefore n G" H 1 / 2 (r). 
We still need to do some work on H ±1//2 (r) related to the normal direction. What follows 
is one of several possible ways. Let Cp := jf| J r x be the center of mass of T and let 
m(x) := x — Cr- Note that m G HV 2 (r), 

/mn = / div(x-cr) = d|fi_|, and / e • m = Ve G P (T). (4.1) 
Jr Ju- Jr 

Proposition 4.1. The decomposition 

H^^r) = Hy 2 (r) © span {m}, with H^r) : = {£ G H 1 / 2 (r) : jf £ ■ n = 0}, 

is stable. Moreover, i"b(r) C Hy 2 (r). 
Proof. Given £ G H 1 / 2 (r), we can decompose 

£ = & + cm, c:=-±— [ £-n = — ^ / £ • n, £o G Hy 2 (r). 

d\tt_\ J r J r m-nJ r 

Since |c| < (rld -1 ^-!" 1 ||£||r, the stability of the decomposition follows readily. The 
inclusion P (T) C H* /2 (r) is a simple consequence of the divergence theorem. □ 
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Proposition 4.2. The decomposition 

ET^r) = H-y 2 (T) ©span{n}, with H^ 2 {T) := {A G H _1 / 2 (r) : (A, m) r = 0}, 

1 /2 

is stable. Moreover, Pq(T) C Hm (r). 
Proof. Given A G H~ 1//2 (T), we decompose 

A = A + cn, c= i^£ = ^_(A,m> r A„eH^(r). 

(n, m)r d|i2_| 

Since \c\ < c? — 1 | — 1 1| m|| i/2,r ||-^||-i/2,r, the stability of the decomposition follows. The 
inclusion P (T) C H„ /2 (r) follows by (JUJ. □ 

Note that Proposition 14.21 is actually the dualization of Proposition 14.11 given the fact 
that 

H" 1/2 (r) = span {m}° and span {n} = Hy 2 (r) . 
We will also consider the following space: 

H 1/2 (r) := {A G H-V^r) : (A, e) r = Ve G P (T)} = P (T)°, (4.2) 
Proposition 4.3. The following decomposition is stable: 

H-V^r) = (H^r) n H 1/2 (r)) © span{n} © P (T), 
Proof. Given A G H _1 / 2 (r), we decompose A = A + cn + e, where 

(A,m) r 1 > 

c := — and e := — T > (A, e,-)re,-, 

(n,m) r |r| J/ 1 

{ei,...,ed} being the canonical basis of lR d . This decomposition is stable and yields 

A GH- 1/2 (r)nHo 1/2 (r). □ 

Proposition 4.4 (Traces of solenoidal fields). The trace operator 7 : V(M d ) := {u G 
W(M d ) : divu = 0} — > Hn 2 (T) is surjective. 

Proof. If u G V(M d ), then 711 G ~H.l/ 2 (T) by the divergence theorem. Let now £ G Hl/ 2 (T). 
We decompose M. d in the bounded set f2_, surrounded by an annular domain Q b with 
boundary T U S and an exterior unbounded domain Q e with boundary S. We then solve 
the following problems 

-Au +Vp_=0 in fL_, -Au 6 + Vp b = in fl b , 

divu_ = in divu;, = in f^, 

711 = £, on T, 7U 6 = £, on T, 

7U5 = 0, on E, 

u , in f2_, 
u := { u bl in Q b , 
in Vt e . 

Then u G W(M d ) and 7U = £. □ 



and define 
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5 Variational stokeslet distributions 



In what follows we will denote: 

O (u, v) := 2u (e(u), e(v)) c , [tu] := 7 "u - 7 + u. 

In this precise moment we start working only in the three dimensional case. Section [6] 
will also be focused in three dimensional space, while we will leave the modifications for 
two dimensions to Section [7J Given A G H _1 / 2 (r), we look for the solution of 

(u A ,p A ) G W(R 3 \ r) x L 2 (R 3 ) (5.1a) 

-2udiv e(u A ) + Vp x = in E 3 \ T, (5.1b) 

divu A = inM 3 \r, (5.1c) 

[ 7 u A ] = 0, (5.1d) 

[t(u A ,p A )] = A. (5.1e) 

Proposition 5.1. Problem (15.1 p is equivalent to 
u A G W(M 3 ),p A G L 2 (M 3 ), 

:u A ,v)-(p A ,divv) K 3 =(A, 7 v) r VvGW(M 3 ), (5.2) 

(divu A ,g) K3 =0 Vg G L 2 (IR 3 ). 

Problem (15 ,2p we// posed. 

Proof. The equivalence of both problems is a simple application of distribution theory. 
Note that div u A = in IR 3 (since u A does not jump across V) but — 2z/div e(u A ) + Vp A = 
only in M 3 \ T. Well posedness follows from the theory of abstract mixed problems [3], [H 
Theorem II. 1.1] as we next show. By Proposition 12.41 div : W(R 3 ) — > L 2 (IR 3 ) is onto. By 
(12. 2p . we have coercivity 

a R 3(u,u) = 2z/||e(u)|| 2 l 3 > i/||Du||| 8 , 

and the right hand side of the inequality is an equivalent norm in W(R 3 ) by Proposition 
12.11 (see also Corollary 12. 2p . so the bilinear form in the diagonal is coercive in the entire 
space. □ 

Definition. To A G H _1 / 2 (r) we associate (u A ,p A ) the solution of (15.11) and then define 
the single layer potential 

S U X := u A , S P A := p\ 

and the associated operators 

VA := 7 ± u A , K'A := {t(u x ,p x )} = |(t + (u A ,p A ) + 1T(u A ,p A )), 

respectively called single layer operator and adjoint double layer operator (it will take a 
while to see why the name of the latter). 
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Proposition 5.2 (Mapping properties). The following operators are bounded: 



(S u , Sp) 
V 



H~ 1/2 (T) W(M 3 ) x L 2 (R 3 ) 

tr 1/2 (r) -> H 1/2 (r), 
H- 1 /2(r) ^ H _1 / 2 (r). 



Proof. It is a direct consequence of the well posedness of (15.21) . the definitions of the 
operators and (13. 6p . □ 

Proposition 5.3 (Jump relations). Let A G H _1 / 2 (r) and let u A := S U A, p\ := S P A. 
Then 

[ 7 u A ]=0, [t(u A ,p A )]=A, and t ± (u A ,p A ) = =f|A + K*A. 

Proof. The jump properties are given by definition of the potential. The normal stresses 
of the potential can then be computed by adding and subtracting the equalities 

i[t(u A) p A )l = \X and lt(u x , Px )}} = K*A. 

□ 

A particular solution. For A := n, it is simple to check that u A = and p\ := ~Xn- 
is the solution of (15. 2p . Therefore S u n = 0, S p n = — Xri-, and by Proposition 13. 1( b) 

Vn = 0, K*n = — |t~(0, 1) = fn. 

The last identity can also be written (K* — |l)n = 0. 

Proposition 5.4 (Symmetry of V). For all A, fj, G H _1//2 (r) ; 

(/x, VA) r = (A, V/i) r , (A, VA) r > 0. (5.3) 

KerV = span{n} and Range V C Hy 2 (r). 
Proof. We associate A h-> (u A ,p A ) and /x >->■ (u^p^). Then 

(/x, VA) r = (^,7 u A)r = a R 3(u M ,u A ) - (p^, divu A ) R3 = a M 3(u M , u A ). 
This proves (15. 3p . Since Vn = (see above), it follows that 

0= (A,Vn) r = (n,VA) r , 

1 /2 

which, by definition, means that VA G H„^(r). 

If VA = 0, then a R 3(u A ,u A ) = 0. By Corollary 12.21 this shows that u A = 0. Going 
to the Stokes equations, this shows that Vp A = in M 3 \ T. Since p\ G L 2 (R 3 ), this 
implies that p x G span{xn_}. Using the fact that A = [t(u A ,p A )J it follows that A = 
c[t(0, Xn-)} = —en for some c G M. □ 
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Proposition 5.5 (Coercivity of V). There exists > such that 

(A,VA) r >C r ||A|| 2 1/2ir VAeH^(r). 

Therefore Range V = Hl/ 2 (T). 

Proof. Recall that by Proposition 14. 2[ the decomposition H 1//2 (T) = Hl/ 2 (T) © span {m} 

— 1/2 

is stable. If A £ H m (r), this implies that 

||A||-i/ 2 ,r = sup 



Using Proposition 14.41 we can define a right-inverse of the trace operator 7^ : nl/ 2 (T) 
V(R 3 ). Then 

l(-V£)r| = |(A,77^)r| = |a R 3(u A ,7 t ^) - (P\, dh^^l 
vU A ,7^)l < fC r ||e(u A )|| R s||£||i /2 ,r, 



and therefore 

||A||_i Ar < z/C r ||s(u A )|| M 3 VA £ H~ 1 / 2 (r), u A := S U A. 

Finally from the proof of Proposition 15.41 it follows that 

2 

(A, VA) r = a R 3(u A ,u A ) > — ^|| A|| 2 1/2r , 

1 /2 

which proves the coercivity of V. Finally, if £ £ H n (Ij, then we can solve the coercive 
problem 

A£H m 1 / 2 (r) (M,VA) r = (M,Or V/x £ H~ 1 / 2 (r) 
and note that by Proposition 14.21 this implies that VA = £. □ 



6 Variational stresslet distributions 

In parallel with Section |5l we restrict our attention to the three dimensional case. Given 
ip £ H 1 / 2 (r) we look for solutions of 



£ W(R 3 \ r) x L 2 (R 3 ) (6.1a) 

-2udiv e(u v ) + Vp v = in R 3 \ T, (6.1b) 

divu^ = inR 3 \r, (6.1c) 

[ 7 uJ = -<p, (6. Id) 

[tK,pJ = 0. (6.1e) 
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Proposition 6.1. Problem ( 16. ip is equivalent to 

u,e w(i 3 \r),^Gi 2 (M 3 ), 

[7uJ + ^ = 0, 

^^(u^,, v) - (p v , div v) M 3 =0 Vve W(M 3 ) 
(divu^,g) M 3\ r =0 Vg G L 2 (IR 3 ) 

Problem (16. 2p we// posed. 



(6.2) 



Proof. Equivalence of (16. ip and (16. 2p is a simple application of distribution theory. Taking 
Lip G W(M 3 ) such that IpyLy] = — ip, we can rewrite the problem in terms of the variable 



u 



L(p e W( 



The resulting problem is associated to the same bilinear form as 



problem (15.21) and is therefore well posed. 



□ 



Definition. To (p G H 1 / 2 (r) we associate (u^p^,) the solution of (16.11) and then define 
the double layer potential 

B u (p := u^, D p (p:=p v 

and the associated operators 

W(p := -t ± (u ¥ ,,p ¥ ,), Kcp := {-fu^} = \ (j + u v + 7~u^) , 
respectively named hypersingular operator and double layer operator. 
Proposition 6.2 (Mapping properties). The following operators are bounded: 



(D„,D P ) 
W 
K 



H 1/2 (r) W(M 3 \ r) x L 2 (M 3 ) 
H 1 / 2 ^) H- 1 / 2 ^), 

H 1/2 (r) -> H 1/2 (r). 



Proof. It is a direct consequence of the well posedness of (16. 2p . the definitions of the 
operators and (13. 6p . □ 

Proposition 6.3 (Jump relations). Lei <p G H 1 / 2 (r) and /ei := T) u tp, p v := T) p (p. 
Then 

[7uJ = [t(u v ,?v)] = 0, 7 ± u^ = ±|v? + Kc^. 



Proof. It is a straightforward consequence of the definitions. 



□ 



A kernel. For c G M, the fields u = — Xn_ c and p = solve (16. ip with (p = c £ Atr 
7 At. This leads to the identities: 



D p c = 0, 



Kc 



2 ' 



Wc = 0, Vc G At r . 



(6.3) 
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Proposition 6.4 (Symmetry of W). For all ip, if) G H 1 / 2 (r), 

(W<^V}r = W,^} r , (W^,^) r >0. (6.4) 

Also 

KerW = At r and RangeWcAtp. 
Proof. Let (u^p^) := (D u cp, D p <£>) and tty := D u t/>. Then 

(Wy>, V)r = (t(u^,^), [7U^]) r = (t~(u v ,p,,), 7~tty) r - (t + (u^,p (/ ,), 7 + u^) r 
= aK3\r(U(p, ity). 

This equality implies symmetry and positive semi-definiteness. 

We already know from (16.31) that Atr C KerW. If Wy? = 0, then, by the previous 
identity, it follows that e(u ¥ ,) = in IR 3 \ V. By Lemma [2.51 it follows that u v ,|n_ G At 
and u|n + = 0. Therefore (p = — [711^] = — 7~u G Air- 

Finally, if cp G H 1 / 2 (r) and c G At, then defining (u^p^) as usual, 

(Wv,c) r = -(t"(u (/ ,,p ¥ ,),7c) r = -2i/(e(u v ), e(c)) n _ + (p^,divc) n _ = 0, 
which proves that Wip G At p. □ 
Proposition 6.5 (Coercivity of W). There exists Cp > snc/i i/iai 

(W^>r>Cr||^| 2 /2 ,r G H^ 2 (r) := : ^ • c = Vc G At r }. 

Therefore Range W = Alp. 

Proof. Proceeding as before, with = Dip, we prove that 

(W<p,<p) T = 2v\\e{u (p ) ' 



Using Proposition 12. 7j Korn's inequality and a compactness argument, it is easy to prove 
that 



lull 2 := llefu 



6 1 r 2 

+ 2_. / [7 U 1 ' °i where At = spanjci, . . . , c 6 }, 



is equivalent to the usual norm in W(M 3 \ T). However, if cp = -[711] G H^ 2 (r), then 
Hull = ||e(u)||ro3\ r and thus 



(Wip,ip) r >C r \\uJl 

Finally 

||^||i/2,r < ||7 + u ¥ ,||i/ 2 ,r + h^Hi/ap < CUu^m^ 

and the proof of the coercivity estimate is finished. To show that the range of W is At r , 
follow the steps of the end of the proof of Proposition I5.5[ using now the decomposition 

r l/2 



H 1 /2(r) = H^ 2 (r)©At r . □ 
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Proposition 6.6. Single and double layer potentials are orthogonal with respect to the 
semi-inner product <%3\ r ( • , • ), that is, 

v = S n A, u = D u (p =^ a K 3\ r (u, v) = 0. 

Proof. Taking v = S U X G W(R 3 ) as test function in problem (I6.2p the result follows. □ 

Proposition 6.7 (K* is the transpose of K). 

(A, K^) r = (K*A, cp) r VA G H" 1 / 2 ^), <p G H J / 2 (r). 

Proof. Let (u x ,p x ) : = (S U A, S p A) and (u^p^,) := (D u (p,D p (p). By Proposition ESI it 
follows that 

(t _ (u A ,p A ),7 _ u (/ ,) r = a n _(u A ,u^) - (p x ,diyu (p ) a _ = a n _(u x ,u v ) 

= -a n+ (u A , u v ) = -a n+ (u A , u^) + (p x , div u^)q + 
= (t + (u A ,p A ),7 + u^) r . 

However, this equality is equivalent 

(t + (u A ,p A ),-[7uJ) r = (|t(u A ,p A )],7~u^) r , 

which, by the jump properties (Propositions [5T31 and 16.3]) . can be equivalently written as 

(_I A + K'A, <p) r = (A, -\<p + K^) r , 

which, once more, is equivalent to the statement. The reader who has followed the proof 
in detail will have realized that this transposition theorem ends up being equivalent to the 
orthogonality property (Proposition I6.6[) . which asserts that energy of single and double 
layer potentials remains separated. □ 

Proposition 6.8 (Representation formula). Let u G W(R 3 \ V) and p G L 2 (R 3 ) satisfy 

-2z/div e(u) + Vp = in R 3 \ T, 
div u = in R 3 \ V. 

Then 

u = S tt [t(u,p)]-D u [7u] and p = S p [t(u,p)] - D p [ T u]. (6.5) 

Proof. Let A := [t(u,p)J and £ := [711]. Define then v := S n A — D, u £ and q = S p A — D p £. 
By definition of the layer potentials, this pair satisfies 

-2z/dive(v) + Vq = in R 3 \ T, 
divv = inR 3 \r, 

[TV] = t 
[t(v,g)] = A. 

However, (u, p) satisfies the same equations. This problem has at most one solution (this 
is part of Proposition 15.11 and 16. ip . Therefore, u = v and p = q, which proves the 
result. □ 
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7 Two dimensional hurdles 



There are two new difficulties in the two dimensional case: 

• The exterior normal stress has zero average on the boundary T (see Proposition 
13.1( d)). This is also reflected in the fact that the fundamental solution does not 
decay and a condition of zero average on the density is needed to ensure decay of 
velocity and pressure at infinity. 

• On the other hand, the associated weighted Sobolev space contains constant func- 
tions -Proposition 12.1} Corollary 12.21 and formula (12. ip -. This fact leaks a two 
dimensional kernel (constant velocity fields with no pressure) in the set of solutions 
in free space. The associated transmission problems that we use to define the layer 
potentials are solvable up to constant velocity fields. Since the integral represen- 
tations of the potentials show decaying velocity fields, we have to choose the right 
one from the variational formulation. The way we are going to do it is by looking 
around with the integral formulation (this is a Deus ex machina moment in the 
theory), picking up an average and imposing it on the solution. There might be 
simpler ways. (For instance, in the Laplace equation, the average around can be 
shown to be zero, which is how the right solution is chosen.) 

As a byproduct of all of this, the representation formula shows that the jumps of the 
Cauchy data represent the solution up to a constant velocity field. This fact ends up 
being the variational form of the Stokes paradox. 



That Fredholm feeling. If we try to define the layer potentials with the transmission 
problems 

(u,p) G W(R 2 \T) x L 2 (R 2 ), (7.1a) 
-2u div e(u) + Vp = in R 2 \ T, (7.1b) 
divu = inR 2 \r, (7.1c) 

bu]=<P, (7- Id) 

[t(u,p)] = A, (7.1e) 

we will be instantly annoyed by the existence of a two dimensional kernel: (u ro ,0) G 
P (JR 2 ) x {0} is, unfortunately, a homogeneous solution of (17. ip . Anyone slightly aquainted 
with Fredholm properties will be ready to see that we need to ground the solution (in 
order to find a unique solution), and that the data (<p, A) G H 1 / 2 (r) x H" 1 / 2 ^) will 
be forced to satisfy a compatibility condition. The second requirement is easy: it falls 
entirely on the demand that A G H (r), which is necessary for existence of solution 
by Proposition 13. 1( c)- (d) . We will ground the solution by integrating around the domain. 
We thus define 

j(u) := / u(x)dS(x), 

where H = dB(0; R) and Q~ C -8(0; R). The value for this surrounding integral will be 
taken from some moments on the two densities. This is where we bring some information 
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from the future (from what it will be when we 'go all integral'). With the matrix valued 
kernels (see ahead in (19.51) ) 



E„(r) := 



— ( logr 1 1 + — r <g> r J , T u (r; n) := — r(r ■ n)(r <g) r), 

4:7iu \ r l J 7rr 4 



we construct two matrix-valued functions on V (generated from S) 
B s (y) := ^ E u (x - y)dH(x), B D (y) := T u (x - y; n(y)) dS(x), y G T, (7.2) 



and then two R 2 -valued linear functionals 

4(A) := (B s , A) r , £ D (tp) := J B D (y)<p(y) dT(y). 

The grounding condition for problem (17.1ft will then be 

j(u)=£ s (A)-£ D M, 

and the space 



{u G W(R 2 ) : j(u) = 0} 



(7.3) 



(7.4) 



(7.5) 



will play the role of the ground zero. From this moment on, the narrative is going to be 
extremely (while not strikingly) similar to that of Sections and O Since at this point 
the reader will be used to these arguments, we are going to be bolder and to define both 
layer potentials at the same time. Most proofs of this section are straightforward copies 
of those in Sections [5] and El 



Proposition 7.1. Problem (17. ip is equivalent to 
uGW(R 2 \r),pGi 2 (R 2 ), 

M = <p, 

a R 2\ r (u,v) - (p,divv) M 2 = (A,7v) r Vv G W(R 2 ) 

(divu,g) M2 =0 Vg G L 2 (R 2 ). 

The modified problem 

uG W(R 2 \r),pG L 2 (R 2 ), 
[7u] = <f, 

j(u)=4(A)-£dM, 

ai;2\ r (u,v) - (p, divv) K 2 = (A,7v) r Vv G W(R 2 ) 

(divu,g) M 2 =0 Vg G L 2 (R 2 ). 



(7.6) 



(7.7) 



is well posed. Therefore, problem (17.61) has a solution that is unique up to elements of 
P (R 2 )x{0}. 
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Proof. The equivalence of (17. ip and (I7.6P is, as usual, a simple exercise in distribution 
theory. Note now that 

W(R 2 ) = W*(R 2 ) © P (R 2 ). 
Since A G H (r), then the decomposition (17.81) shows that 



(7i 



a R 2(u A ,v) - (p A ,divv) R2 = (A,7v) r Vv G W*(R 2 ) 

if and only if 

a R 2(u A ,v) - (p A ,divv) R2 = (A,7v) r Vv G W(R 2 ). 
We next get rid of the non-homogeneous side conditions. We start constructing 

u;GW(R d \T), 7 X = <p, u;^0infi + , 
so that f/yu°] = (f and j(u° ) = 0. We then propose 

u = u* + u; + ^(4(A)-£ D M), 

where 

u* G W„(R 2 ),p G L 2 (R 2 ), 

a M2 (u*,v) - (p,divv) R2 = (A, 7v) r - a n _(u° , v) Vv G W„(R 2 ), 
(divu*,g) R2 =0 Vg G L 2 (R 2 ), 

as the solution of (17. 7p . That problem (17.91) is well posed follows from the theory of mixed 
problems and the arguments we next give. By Corollary 12. 2[ ||e(-)l|iR a ^ s equivalent to 
the weighted Sobolev norm in W Vc (R 2 ), which means that the diagonal bilinear form a^2 
is coercive m W*(R 2 ). By Proposition El div : W(R 2 ) L 2 (R 2 ) is onto and therefore, 
using (J72D, it follows that div : W*(R 2 ) ->■ L 2 (R 2 ) is also surjective. □ 



(7.9) 



Two (or four) potentials, four operators. Problem (17.71) defines a bounded linear 
map HV2(r) x H 1/2 (r) W(R 2 \ T) x L 2 (R 2 ). We are then allowed to see this as a 
matrix of operators 

A 



u 




S u 


-D u " 




A 


. p - 




S p 









[ 



] 



where the second expression looks at (u,p) as a joint entity, and not as separate fields. 
The integral operators are defined by taking averages 

[V K] :={{ 7 -}}[S u D n ], [K* W ] := ft( • )| [ S -D ] . 

Proposition 7.2 (Mapping properties). The following operators are bounded: 



(S n 


S p ) 


Hq 1/2 (r) ->■ W(R 2 ) x L 2 (R 2 ), 


(D„, 


D P ) 


H 1/2 (r) W(R 2 \ r) x L 2 (R 2 




V 






K* 






K 


H 1/2 (r) -> H 1/2 (r), 




W 





19 



Proposition 7.3 (Jump relations). For all A G H 1 ^ 2 (r) and ip G H 1 / 2 (r) ; 

[ 7 S U A] = 0, [yD u ip] = -if, 
[t(S u A, S P A)] = A, [t(D ttV >, D p¥ »)] = 0, 
t±(S M A, S P A) = T ±A + K*A, t ± (D u ^, D p cp) = ±\<p + K^>. 

Simple solutions. With A = n and ip = 0, we get 

S u n = 0, S p n=- X n-, Vn = 0, -§n + K*n = 0. (7.10) 
With A = and <p = c G Ai r = 7.M, we get 

D»c = -xn_c, D p c = 0, Wc = 0, ±c + Kc = 0. (7.11) 

— 1/2 

Proposition 7.4 (Properties of V). For all A,/iG H (T) ; 

( M , VA) r = (A, VA*) r , (A, VA) r > (7.12) 
and there exists Cr > such that 

(A, VA) r > C r ||A||! 1Ar VA G H m 1 ^ 2 (r) n H 1/2 (r). 

Also 

KerV = span {n} and Range V C H^ 2 (r). 

Proof. The proof of Propositions 15.41 and 15.51 (same results in the three dimensional case) 
are valid. □ 

The range of V. Identifying the range of V requires some additional work, given the 
fact that we have been defining V on H (r), as opposed to having it defined on the 
entire H _1//2 (r). To do that we consider the solution of the coercive problems (J = 1,2) 

^ eq G H-V2(r), «,e,) r = l, 
(^ eq ,V^) r = V/i G H -1 / 2 (r), 

where {e!,e 2 } is the canonical basis of K 2 . These densities are a basis of the space of 
equilibrium distributions for the Stokes problem (see [7] and [6]). They are the Stokes 
equivalent of the electrostatic equilibrium distribution that leads to the definition of log- 
arithmic capacity. Using a coercivity argument (see the proof of Proposition 15. 5p . it is 
easy to see that 

V : H 1/2 (r) n H m 1 / 2 (r) — ► span {n, <^ q , 2 q }° 

is invertible. 
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Proposition 7.5 (Properties of W). For all (p, ijj G H 1 / 2 (r), 

(W^,V)r = (Wif>,<p)r, (W<p,tp) r >0 (7.13) 
and there exists Cr swc/i t/icrf 

(W^¥»)r>C7r||v||?/2,r V^> G H^ 2 (r) := {| : ^-c = Vc G M r }- 

KerW = .Mr and RangeW = A^. 

Proof. We follow the proof of Proposition 16.41 Let (u^,^) := (D u (p,D p cp) and : = 
T) u if). Then we can prove that 

{W<p, il>) r = a R 2\ r (u ¥ ,, u^), (7.14) 

which implies symmetry and positive semi-definiteness (I7.13p . 

Formula (l7TTTj> shows that Mr C Ker W. If W<£ = 0, then, by fl7TT4"j) . it follows that 
e(u^) = in M 2 \ T. By Lemmas 12.51 and 12.61 it follows that u v |n G M and u ¥ ,|q + G 
Pq(Q+). Therefore ip = — [711^] G Mr- The remainder of the proof of Proposition 16.41 
can be applied verbatim. 

To prove coercivity, we proceed as in the proof of Proposition 16.51 We first need to 
prove that 



lull 2 



3 1 f 2 

(u)|| 2 l 2\ r + |j(u)+£ D ([7u])| 2 + ^ / [7u]-q , where M = span{c x , c 2 , c 3 }, 

1=1 1 Jr 



is a norm in W(IR 2 \ T). This can be done using Lemmas 12.51 and 12.61 Using Korn's 
inequality and a compactness argument, we can show that this is an equivalent norm in 
W(IR 2 \ T). However, |||u||| = ^(u^H^p if <p G H^ 2 (r). The remainder of the proof of 
Proposition 16.51 can be then applied to show coercivity. 

Finally, using coercivity and the fact that RangeW C M r , it is easy to see that both 
sets are equal. □ 

Proposition 7.6 (K* is the transpose of K). 

(A, K^) r = (K*A, <p) r VA G H^r), <p G H 1 / 2 (r). 

Proof. We first prove that single and double layer potentials are orthogonal with respect 
to the semi-inner product a R 2\ r ( • , • ), that is, 

u = S U A, v = D u <p =^ a K 2\ r (u, v) = 0. 

This is done by reversing the roles of u and v in (j7.7j) . given the fact that S U A G W(M 2 ) 
can be used as a test function. The remainder of the argument follows the proof of 
Proposition 16.71 □ 
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Proposition 7.7 (Representation formula). Let u E W(R 2 \ T) and p G L 2 (R 2 ) satisfy 

-2z/div e(u) + Vp = in R 2 \ T, 
div u = in M 2 \ T. 

Then there exists Uoo 6 Po(^ 2 ) such that 

u = u 00 +S tt [t(u,p)]-D w [ 7 u] and p = S p [t(u,p)] - D p [ 7 u]. (7.15) 

Proof. See the proof of Proposition 16 .81 The argument uses the fact that the unique solu- 
tions to (17.11) with vanishing right hand side are the elements of Po(^ 2 ) x {0} (Proposition 

ED. □ 



8 Review #2: Laplace and Lame potentials 

In this section we give a fast review of some known properties of the layer potentials for 
the Laplace and Lame (or Navier, or homogeneous isotropic linear elasticity) equations. 
The variational approach (with weighted Sobolev spaces, as in Sections El El and [7]) can 
be found in [T3] , although the fine print of relating integral forms (what we will do in the 
remainder of the paper) is missing. A purely integral theory is developed in [5] and [121 
Chapter 8 & Chapter 10]. Note that the variational theory gives fast and elegant proofs 
of mapping properties in the basic variational setting but needs some additional work in 
order to show that the variational definitions correspond to weak integral potentials. 



Laplace potentials. For questions of easiness of reference, we will make copies of the 
layer potentials for the Laplacian so that they act componentwise on vector-valued den- 
sities. The single layer and double layer potentials for the Laplacian are defined for 
xGt^Tby 

(S A A)(x) := (E A (x--),A)r, (D A V )(x) := J T A (x - y; n(y)) V (y)dT(y), 
where 

T A (r; n) := — - i ~(r ■ n)I, 

contain diagonal copies of the Laplacian monopole and dipole distributions, and r = |r|. 
For some mapping properties, we will use the spaces H 1 1 oc (R 2 ). These spaces are not 
normed spaces but are easily seen to be metrizable, by using a sequence of cut-off functions. 

Proposition 8.1. For all (A, ip) e ET 1 / 2 ^) x H 1 / 2 (r) ; the function u := S A A - D A ip 

is a solution of 

-Au = mR d \T, (8.1a) 
[ 7 u] = <p, (8.1b) 
[(Du)n] = A. (8.1c) 
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Moreover, the following mappings are continuous: 

( H~ 1 / 2 (r) — > W(K 3 ), d = 3, 
D A : H 1/2 (r) -)• W(R d \T), S A J H 1/2 (r) ->■ W(M 2 ), d = 2, (8.2) 

I H~ 1/2 (T) ->■ Hj 0C (R 2 ), d = 2. 

Lame potentials. For reasons that will become apparent at the time of comparing 
potentials with the Stokes problem, all Lame potentials will be written in terms of one of 
the Lame parameters \x and the non-physical quantity: 

A, " 



A + 2/i 

The layer potentials are defined for x e ~R d \ V by 

(S L A)(x) := (E L (x--),A)r, (D L y>)(x) := J T L (x - y; n(y)) V (y)dT(y), 
where (see pH (2.2.21)] for the definition of T L ) 

^« : =4^i5i(( 1 + A ^ I+ ( 1 - A ^ r ® r )- 
TL(r; n) := 2]r (d - 1) (^ ((r ' n)I " r ® n + n 8 r ) + (1 " A) ^ (r ' n)r ® r ) • 

Note that the definition of the Lame dipoles T L ( • ; n) can be found with the process 

l d 9d,nMT L (.;n) T d = -<r L (E L d)n = - (2//e(E L d)n + A div (E L d) n 

= -^(2Ae(E L d)n+ (1 - 2 A)div (E L d) n) , 

where d acts as a dipole direction and n as the stress direction. 

Proposition 8.2. For all (A, ip) G H" 1 / 2 ^) x H 1 / 2 (r) ; the function u := S L A - D L <p is 

a solution of 

-div cr L (u) = m E^T, (8.3a) 
[7u] = <p, (8.3b) 
[<r L (u)n] = A. (8.3c) 

Moreover, the following mappings are continuous: 

( U-^ 2 (T) -> W(M 3 ), d = 3, 
D L : H 1/2 (r) W(M d \T), S L : I H 1/2 (r) -> W(M 2 ), d = 2, (8.4) 

{ H-V^D^HL^ 2 ), d = 2. 

Note how, if we write cr A (u) := Du, Propositions I8.ll and 18.21 are almost identical 
twins. 
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9 Stokes from Laplace and Lame 

To close this exposition of the Stokes layer potentials, we show that their integral expres- 
sions give the same functions as those defined in Sections M and UJ We will avoid having 
a redundant name for two entities that end up being the same, by ignoring that we ever 
defined the variational potentials in the previous sections. Consider first the pressure part 
of the potentials (as usual, x6l d \r) 

(SpA)(x) := (A,e p (x--))r, (D pV )(x) := J t„(x - y; n(y)) • ^(y)dr(y), (9.1) 



where 



1 

2n(d~l)r d ~ 2n(d-V 



tp(r;n) := 2^T) 2z/ (^ r0r- ^) n - 

It is easy to see that both potentials define functions with C°° (M d \ T) components. 
Proposition 9.1. The following identities hold in R d \ T for arbitrary A e H _1 ' 2 (r) and 

S P A = -div S A A, D p ip = -(2i/)div D A (p. (9.2) 
Therefore, the following maps are continuous: 

( H- x / 2 (r) -> L 2 (R 3 ), d=3, 
Dp : H 1 ^ 2 (r) — >■ L 2 (IR d ), Sp : i Ho 1/2 (r) ^ L 2 (M 2 ), d = 2, (9.3) 

1 H-V2(r)^LL(]R 2 ), d = 2. 

Proof. The first part follows by a simple computation. Mapping properties are a simple 
consequence of Proposition 18. II and of f)9.2p . □ 

The velocity part of the Stokes potentials is defined by respective superposition of 
stokeslets and stresslets: 



(S„A)(x) := (E u (x--),A) r , (D u¥ >)(x) := / T u (x - y; n(y)) V (y)dr(y), (9.1) 
where 



r 



E u (r) := (V(r)I + i-r ® r) , (9.5a) 



4vr(d- 

Tu(r; n) := 2n(d- 1) r^+ 2 (r ' n) r ® r ' (9 ' 5b) 
For some forthcoming arguments, it will be useful to have the rotlet distributions at hand 

(R^)(x) := f M(x - y; n(y)) A(y)dr(y), (9.6) 
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r 



where 



M ( r ' n ) := ^T7~a — Tw (r ® n - n ® r). 

lix{d — l)r a 



Moreover 



Proof. Note first that 



A simple computation shows that in R d \ T: 

div Rep = -div D A tp = (2v)- 1 T) p ip G H 1/2 (r). (9.7) 

Proposition 9.2. T7ie following mappings are continuous: 

( H _1 / 2 (r) — >■ W(IR 3 ), d=3, 
D u : H 1 / 2 ^) ->■ W(lR d \ T), S u : H " 1/2 (r) ^ W(i 2 ), d=2, (9.8) 

I H-nr) ^ HLOn d = 2. 

[ 7 D m¥ >] = -y> G H 1 / 2 (r). (9.9) 

S L = /i- 1 AS A + (l-A)z//i- 1 S u , (9.10) 

which transmits the mapping properties of the Laplace and Lame single layer potentials 
(Propositions 18.11 and 18. 2p to S u . Note that this includes the fact that [7S U A] = 0. We 
next notice that 

D L = A D A - A R + (1 - A)D U = D u + A (D A - R - D tt ). (9.11) 

Since D L is bounded H 1 / 2 (r) — > W(R d \T) for all mathematically valid choices of the Lame 
parameters (and then at least for all < A < 1), it is clear that D u and D A — R — D n have 
the same mapping property. Note that this implies that the rotlet operator is bounded 
between the same spaces as well. Taking the jump of the trace on both sides of ( 19.1 ip 
and using Propositions 18.11 and 18.21 it follows that 

_^ = _ A(p - A^Rcpj + (1 - A)[iD u <p], 

again for all < A < 1, and therefore ( 19. 9 p holds. Note that, in passing, we have proved 
that [7Ry>] = 0, and therefore R : H 1 / 2 (r) -»■ W(R d ) is bounded. □ 

Proposition 9.3. Let (A, ip) G H" 1 / 2 ^) x H 1 / 2 (r) and let 

u := S M A - B u (p, p : = S P A - D p (p. 

Then 

(u,p) G Ul oc (R d \ r) x Ll c (R d ), (9.12a) 
-2z/div e(u) + Vp = in R d \ T, (9.12b) 
divu = mR d \T } (9.12c) 
[ 7 u] = ip, (9.12d) 
[t(u,p)] = X. (9.12e) 

Moreover, (u,p) G W(M d \ T) x L(M d ) if d = 3, or if A G H 1/2 (r) and tZ = 2. Fma%, 
w/ien d = 2, condition (17.41) zs satisfied. 
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Proof. Propositions 19.11 and 19.21 prove (19.12aj) and (19.12d|) . The differential equations 
(19.12bj) and (19.12dj) can be proved to hold by a direct (but cumbersome) computation. 



Let now A G H 1 ^ 2 (r), take fi = v and A = in the Lame equations (therefore A = |), 
define 

u := S U A, p := S P A, u L := S L A, u A := S A A, 
and note that by (19.101) and Proposition 19.11 

u = 2u L — z/ _1 u A and p = — divu A . 



By Propositions 18.11 and 18.21 and the definitions of the normal stresses, it follows that 

2z/(£(u L ),Dv) Kd = (A, 7 v)r, (Du A ,Dv) Rd = (A, 7 v) r , Vv G Hl omp (R d ), 

where the subscript 'comp' is used to denote compact support. Therefore 

([t(u,p)],7v)r = 2i/(e(u),Dv) R d - (p,divv) R d 

= 4i/(e(u L ), Dv) Rd - 2(6(u A ), Dv) Rd + (div u A , div v) Rd 
= (A,7v) r - (Du A , (Dv) T ) Rd + (divu A , div v) Rd , 

for all v G H* omp (R d ). However, it can easily be seen with a density argument and 
differentiation in the sense of distributions that 

(Du A , (Dv) T ) Rd = (divu A ,divv) Rd Vu G Hl c (R d ), v G H^K*). 

This proves that 

[t(S u A,S p A)] = A VAGH- 1 / 2 ^). 
We start afresh with <p G H 1 / 2 (r), take general A and fi, and define 

u:=D u ^>, p:=D p <p, u L := D L cp, u A := DV u R := R^, 

and note that by (19. 7p and (19. lip , it follows that 

div u R - div u A = is' 1 p and u L = A u A - A u R + (1 - A)u. (9.13) 

Noticing that 

Afi^cr^iy) = 2Ae{y) + (1 - 2A)(divv)I, 
using (19.131) and the fact that div u = 0, we can expand 

A /U -V L (u L ) = 2A 2 (e(u A ) - e(u R ) - e(u) - (divu A )I + (div u R )I + (div u)l) 
+A (2e(u) + (divu A )I- (divu R )I- (3divu)l) + (divu)I 
= 2A 2 (e(u A )-£(u R )-s(u) + z/- 1 j9i) +A(2e(u)-u- 1 pI). 

Taking the jump of the normal components on both sides of the previous formula, it 
follows that 

= /x-VW] = 2A[(6(u A - u R - u) + ^I)n] + iz-^tCu.p)]. 
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Since this results holds for all < A < 1, then [t(u,p)J = 0. 

If we now integrate u on the curve/surface H = dB(0;R) and note that -with the 
notation of (Q and (Pl) - 




-y;n(y))^(y)dr(y))dS(x) 
-y;n(y))dS(x))cp(y)dr(y) 



= (B s , A) r - y B D (y)^(y)dr(y) = £ S (X) - £ D (<p), 
then we get condition (I7.4p in the two dimensional case. □ 



In conclusion. Proposition 19.31 shows that the potentials defined with integral formu- 
las (19.11) and (19. 4p are the same as the potentials defined through variational problems 
(Sections El [6] and [7]). Furthermore, formulas (19.21) . (I9.10p . and (19.111) give simple expres- 
sions of the Stokes potentials in terms of Lame and Laplace potentials, with the fleeting 
presence of the rotlet (19.61) . Our story ends here, but this is far from all that can be done 
with the theory of these operators. The formulas (19.101) and (19. lip (and the parametric 
analysis in terms of the non-physical quantity 0) give a fast track transfer of many results 
from the better understood theory of Lame and Laplace layer potentials to the Stokes 
world. 
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